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Abstract. Ferrario & Wickramasinghe (2006) explored the hypothesis that the magnetic fields of neutron stars are of fossil 
origin. In this context, they predicted the field distribution of the progenitor OB stars, finding that 5 per cent of main sequence 
massive stars should have fields in excess of 1 kG. We have carried out sensitive ESPaDOnS spectropolarimetric observations 
to search for direct evidence of such fields in all massive B- and O-type stars in the Orion Nebula Cluster star-forming region. 
We have detected unambiguous Stokes V Zeeman signatures in spectra of three out of the eight stars observed (38%). Using a 
new state-of-the-art Bayesian analysis, we infer the presence of strong (kG), organised magnetic fields in their photospheres. 
For the remaining five stars, we constrain any dipolar fields in the photosphere to be weaker than about 200 G. Statistically, 
the chance of finding three ~ 1 kG fields in a sample of eight OB stars is quite low (less than 1%) if the predictions of Ferrario 
& Wickramasinghe are correct. This implies that either the magnetic fields of neutron stars are not of fossil origin, that the 
flux-evolution model of Ferrario & Wickramasinghe is incomplete, or that the ONC has unusual magnetic properties. We are 
undertaking a study of other young star clusters, in order to better explore these possibilities. 
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INTRODUCTION 

The origin of the magnetic fields of massive stars, present 
at all evolutionary stages (PMS, MS, post-MS [1-3]), 
remains an open question. Two possible models may 
explain the presence of magnetic fields: 

(i) In the dynamo model, the field is generated by a 
dynamo mechanism, occurring classically in convective 
regions or induced by strong shear during differential 
rotation. 

(ii) In the fossil model, the field is a remnant from 
a dynamo active during a previous evolutionary phase, 
or swept up from the interstellar medium (ISM) during 
star formation. This scenario implies that the field must 
somehow survive the various structural changes encoun- 
tered during stellar evolution. The magnetic flux is usu- 
ally assumed to be conserved to some extent. 

The dynamo model reproduces well the characteristics 
of late-type main sequence stars and giants. However, it 
fails to explain the fields of magnetic Ap stars, as their 
envelopes are entirely radiative. Some models of convec- 
tion in the small convective core of those stars have been 
put forward, but they still have fundamental difficuties 
reproducing the observed fields [4]. Their simple mag- 
netic geometries, lack of significant mass-field strength 
or period-field strength relation, and the fact that the ob- 
served characteristics of Herbig star magnetic fields [e.g. 
5-9] are qualitatively identical to those of Ap stars, point 



toward a fossil origin. 

Furthermore, the incidence, geometries and strengths 
of white dwarf magnetic fields are at least qualitatively 
compatible with evolution from Ap-Bp stars, suggesting 
that the fields of white dwarfs may also be of fossil origin 
[10]. 

Similar questions can be asked about neutron stars. 
Neutron star field strengths, inferred from spin down 
rates of radio pulsars, are in the range of 10 n -10 14 G 
A particular class of pulsars called the magnetars, which 
include the anomalous X-ray pulsars and soft gamma 
repeaters, have fields that range up to 10 14 -10 15 G. The 
two possible scenarios that have been put forward to 
explain the existence of these fields are again the fossil 
and dynamo mechanisms. 

The fossil hypothesis implies that neutron star fields 
come from the progenitor OB star fields which have 
survived the post main sequence and the core collapse 
phases. There is some observational evidence that neu- 
tron stars may evolve from stars as massive as 45M (;) 
[11, 12]. The only two known O stars with directly 
measured magnetic fields are 0'OriC and HD 191612 
(~ 40 M ), with dipolar field strengths ranging between 
1000-1500G [13, 14]. The magnetic flux of these stars 
(~ 10 27 Gcm 2 ) is of the same scale as the magnetic flux 
of the highest field magnetars. 

Provided that OB star fields are remnants from the 
ISM, the fossil hypothesis could provide a powerful ex- 



planation of the wide range of magnetic fields present 
in neutron stars, and may also explain the super strong 
fields seen in magnetars. 

In this fossil model, the properties of neutron star mag- 
netic fields are a function of the field properties of the 
progenitor OB stars, plus any physics that occurs during 
post main sequence evolution to alter those fields. De- 
pending on the relative importance of these two ingredi- 
ents, the fields of neutron stars might be nearly identical 
to those of OB stars (but stronger, of course), or rather 
different. 

On the other hand, it has been suggested that neutron 
star magnetic fields could instead be generated by a dy- 
namo mechanism taking place during the core collapse 
itself, and induced by differential rotation. Present stud- 
ies [15] assume that any primordial fields present in the 
progenitor star are weak enough to be expelled by the 
dynamo process. However, if the initial field is strong 
enough, the evolution will be different, as this field is 
likely to interfere with the differential rotation and there- 
fore with the dynamo process itself [16]. 

Hence, in both cases, and the likely combination of 
them, a primordial field present during the formation of 
a neutron star will play a fundamental role. Knowing the 
magnetic properties of the progenitor OB stars would 
therefore be an important asset for constraining models 
of stellar evolution leading to neutron star birth. 



OBSERVATIONS 

Magnetic fields can be directly detected in stellar atmo- 
spheres by the means of the Zeeman effect. If the field 
is strong enough, and the spectral lines narrow enough, 
one can directly see the Zeeman splitting of the lines in 
the intensity spectrum. However, if the field is weaker, 
and the lines broadened either intrinsically or by fast stel- 
lar rotation, the splitting is much more difficult to detect, 
even at high spectral resolution. 

In that case, the most effective way to detect the Zee- 
man effect is by looking at circular polarisation signa- 
tures across photospheric spectral lines. Recently, a tech- 
nique called "Least Squares Deconvolution" (LSD) has 
been developed by Donati et al. [17] for extracting a 
mean Stokes V profile, using all the lines present in a 
spectrum simultaneously. This allows for better detection 
limits as it substantially increases the signal-to-noise ra- 
tio. This technique has been proven to be very useful for 
detecting magnetic fields in Ap-Bp stars [18] and late- 
type stars [19]. 

However, magnetic fields in hotter OB stars remain 
difficult to detect. The few photospheric lines present 
in the optical spectrum and the large intrinsic width of 
the lines, worsened by the usual fast rotation of those 
stars, require large-bandwidth and high signal-to-noise 
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FIGURE 1. X-ray efficiency as a function of effective tem- 
perature. The detected stars are circled. Filled symbols are for 
stars with indirect indications of the presence of a magnetic 
field and gray symbols are for confirmed or suspected bina- 
ries. Plotting symbols indicate the following properties: circles 
are for T-Tauri type emission, triangles are chemically peculiar 
(CP) stars, and the diamond star was not observed. The dotted 
line indicates the typical efficiency for massive stars. 



ratio observations to start with, even using LSD. 

For those reasons, the magnetic characteristics of 
the population of neutron star progenitors are currently 
poorly known, and mostly only by indirect ways, such 
as non-thermal radio and X-ray emission and variabil- 
ity, peculiar abundances or cyclical variations of pho- 
tospheric and wind spectral lines (see [19, 20] for a 
review). Nevertheless, those indicators are often in- 
terpreted to show that magnetism may be widespread 
among massive OB stars. 

The advent of a new generation of spectropolarimeters 
such as ESPaDOnS at the Canada-France-Hawaii Tele- 
scope (CFHT) and its twin NARVAL at the Telescope 
Bernard-Lyot (TBL) now allows new investigations of 
magnetic fields in massive OB stars. ESPaDOnS con- 
sists of a polarimetric module located at the Cassegrin 
focus of the CFHT, linked by optical fibres to the high- 
resolution echelle spectrometer. A resolution of 65,000 
for a spectral range covering from 360nm to 1 jim can 
be achieved in a single observation. 

A complete circular polarisation observation consists 
of series of 4 sub-exposures between which the polarime- 
ter quarter-wave plate is rotated back and forth between 
position angles. This procedure results in exchanging or- 
thogonally polarised beams throughout the entire instru- 
ment, which makes it possible to reduce systematic er- 
rors. 
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FIGURE 2. Least Squares Deconvolved profiles for Par 1772 (B2 V) and 1 OriD (B0.5 V). The curves are the mean Stokes I 
profiles (bottom), the mean Stokes V profiles (top) and the N diagnostic null profiles (middle), in solid line for January 2006 and 
dashed line for March 2007. 



The Orion Nebular Cluster 

The first goal of this project was to explore the con- 
nection between magnetic fields and X-ray production in 
massive OB stars. Stellar magnetic fields are well known 
to produce X-rays in late-type convection stars like the 
Sun. However, X-ray emission coming from OB stars 
is often explained by radiative instabilities resulting in 
a multitude of shocks in their winds [21, 22]. 

The Chandra Orion Ultradeep Project (COUP) was 
dedicated to observe the Orion Nebula Cluster (ONC) in 
X-rays. The OBA sample (20 stars) was studied with the 
goal of disentangling the respective roles of winds and 
magnetic fields in producing X-rays [23]. The produc- 
tion of X-rays by radiative shocks should be the dom- 
inant mechanism for the subsample of 9 OB stars with 
strong winds. However, aside from two of those stars, 
all targets showed X-ray flux intensity and/or variability 
which were inconsistent with the small shock model pre- 
dictions (Figure 1). 

For these reasons, we started our investigation with the 
9 OB stars of this young star cluster. They range from 
B3 V (~ 8M ) to 07 V (~ 40 M ), approximately the 
mass range from which neutron stars are thought to be 
formed. 

We conducted spectropolarimetric observations of 8 of 
those massive stars in January 2006 and March 2007 at 
CFHT. Using the LSD technique, we found clear Stokes 
V signatures for 3 stars: the previously-detected 9 l OriC, 
as well as Par 1772 (shown in Figure 2, along with the 
non-detection case 9 1 Ori D), and NU Ori. 



Magnetic analysis 

In order to extract the surface field characteristics from 
the observed Stokes V profiles, we compared them with 
theoretical profiles for a large grid of dipolar magnetic 
configurations, calculated with the polarised LTE radia- 
tive transfer code Zeeman2 [24, 25]. 

We sampled the 4-dimentional parameter space (i, 
/3, <j>, B) which describes a centered dipolar magnetic 
configuration. In such a model, ;' is the projected angle 
of the rotation axis to the line of sight, j3 is the angle 
between the magnetic axis and the rotation axis, is the 
rotational phase and B is the polar field strength. 

For each configuration, we calculated the reduced % 2 
of the model fit to the observed mean Stokes V profiles. 
Assuming that only the phase may change between two 
observations of a given star, the goodness-of-fit of a 
given (z, j3, B)-configuration is expressed in terms of 
Bayesian probability density. This ensures that a good 
magnetic (;, j3, B)-configuration possesses phases that fit 
both observations, as the rotational period is not known 
with enough accuracy to determine a priori the phase 
difference. 

We can determine the probability density of the field 
strength, by marginalising over inclination (;) and obli- 
quity (j3). Then, we extract a 90% credible region for the 
field strength of each star (Figure 3) with the technique 
described by Gregory [26]. 

Our first results show that any dipolar fields present 
in the 5 undetected stars are weaker than ~ 200 G with 
a 90% confidence. The field strengths of the 3 detected 
stars are approximately 1-3 kG. 



0.0025 



0.0000 



Par 1772 
Marginalized posterior probability 



6' OriD 

Marginalized posterior probability 




1000 



2000 3000 
B, (G) 



4000 



5000 



100 



200 300 
B, (G) 



400 



500 



FIGURE 3. Marginalized posterior probability densities for Par 1772 (one of the 3 detected OB stars) and 1 OriD (one of the 5 
undetected stars). The magnetic field strength 90% credible region (filled) is [900, 2900] G for Par 1772 and [0, 160] G for 1 OriD. 
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FIGURE 4. Predicted magnetic field repartition of 8 
randomly-drawn stars according to the predicted distribution 
of Ferrario and Wickramasinghe [27] along with the observed 
distribution of field strengths in the ONC. 



DISCUSSION AND CONCLUSION 

As an illustrative example, we can compare our new ob- 
servational results with the predictions made by Ferrario 
and Wickramasinghe [27] about the magnetic field distri- 
bution of massive stars (8-45 M©) on the main sequence. 

They parametrized the magnetic flux distribution on 
the main sequence x(4>) as the sum of two Gaussians, 
along with the birth spin period of neutron stars. Assum- 
ing a complete conservation of magnetic flux, they have 
calculated the expected properties of isolated radio pul- 
sars. They used the 1374-MHz Parkes Multi-Beam Sur- 
vey of isolated radio pulsars in order to constrain the 
model parameters. They obtained a continuous magnetic 
field distribution in the progenitor OB stars peaking at 



~ 46 G with 5 per cent 1 of the stars having a field in ex- 
cess of 1 kG. 

Of course, our sample contains only 8 stars, but we 
can still make some rough comparaisons. Taking the pre- 
dicted field strength distribution, we assume that it is the 
true parent distribution from which we draw a random 
sample of 8 stars. We define three possible outcomes: [0- 
200] G, [200-1 000] G and over 1000G, with respective 
probabilities derived from the parent theoretical distri- 
bution (Figure 4). According to the multinomial distri- 
bution, the probability of observing the distribution of 
magnetic field strengths observed in the ONC is below 
1%. 

This result might be interpreted, at first glance, to 
suggest that the fields of neutron stars are not of fossil 
origin. However, some points are important to consider: 

1. The sample of stars may not be representative of 
a general parent distribution, as the stars all come 
from the same cluster. This region could be unusu- 
ally magnetic, especially if the fields of the OB stars 
themselves are also of fossil origin from the ISM. 

2. There may be a fossil component to the magnetic 
field origin, but the assumed parent distribution is 
not in fact the true parent distribution because some 
assumptions are incorrect, or some elements are 
missing from the model. Examples of such missing 
physics might be partial flux conservation or the 
influence of dynamo processes during core collapse. 

In order to better explore these possibilities, a larger 
sample of OB stars, from clusters and from the field, 



1 Although the paper states 8%, recalculation based on the detailed 
model distribution provided by L. Ferrario gives 5%. 



must be studied in order to increase the population of 
neutron star progenitors with known magnetic properties. 
Our team has undertaken an extensive spectropolarimet- 
ric study of massive stars in other young star clusters to 
provide these important data. 
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